Discovery of genomic regions associated with grain yield and agronomic traits in Bi-parental populations of maize (Zea mays. L) Under optimum and low nitrogen conditions

Low soil nitrogen levels, compounded by the high costs associated with nitrogen supplementation through fertilizers, significantly contribute to food insecurity, malnutrition, and rural poverty in maize-dependent smallholder communities of sub-Saharan Africa (SSA). The discovery of genomic regions associated with low nitrogen tolerance in maize can enhance selection efficiency and facilitate the development of improved varieties. To elucidate the genetic architecture of grain yield (GY) and its associated traits (anthesis-silking interval (ASI), anthesis date (AD), plant height (PH), ear position (EPO), and ear height (EH)) under different soil nitrogen regimes, four F3 maize populations were evaluated in Kenya and Zimbabwe. GY and all the traits evaluated showed significant genotypic variance and moderate heritability under both optimum and low nitrogen stress conditions. A total of 91 quantitative trait loci (QTL) related to GY (11) and other secondary traits (AD (26), PH (19), EH (24), EPO (7) and ASI (4)) were detected. Under low soil nitrogen conditions, PH and ASI had the highest number of QTLs. Furthermore, some common QTLs were identified between secondary traits under both nitrogen regimes. These QTLs are of significant value for further validation and possible rapid introgression into maize populations using marker-assisted selection. Identification of many QTL with minor effects indicates genomic selection (GS) is more appropriate for their improvement. Genomic prediction within each population revealed low to moderately high accuracy under optimum and low soil N stress management. However, the accuracies were higher for GY, PH and EH under optimum compared to low soil N stress. Our findings indicate that genetic gain can be improved in maize breeding for low N stress tolerance by using GS.


Introduction
Maize (Zea mays L.), a staple crop in sub-Saharan Africa (SSA), is projected to increase in demand (Ekpa et al., 2018), necessitating improved breeding systems to maximize productivity.Over the years, maize productivity in SSA has been limited by the adverse effects of parasitic weeds (e.g., Striga sp.(Kanampiu et al., 2018;Yacoubou et al., 2021)), insect infestations (e.g., fall armyworm (Baudron et al., 2019;De Groote et al., 2020)), and disease incidences (e.g., Maize lethal necrosis (Boddupalli et al., 2020)).Most notably, inadequate soil fertility (particularly low nitrogen) is a primary cause of low maize yields in the region's smallholder farming communities.Modelling studies have established that climate change has the potential to worsen these current maize production constraints (Tesfaye et al., 2015), especially under low soil nitrogen conditions (Falconnier et al., 2020).In this regard, maize breeding holds a significant potential to provide sustainable solutions to the prevailing and projected biotic and abiotic constraints (Eriksson et al., 2018;Ndlovu, 2020).
As nitrogen use efficiency (NUE) is critical for sustainable productivity, there is major interest in developing varieties that perform better in low soil nitrogen stress conditions (Ribaut et al., 2007).However, to harness plant breeding for such stress conditions, an understanding of the performance of maize genotypes under low nitrogen conditions is required.Assessing the performance of maize genotypes under different nitrogen regimes is critical for identifying promising parental lines.Phenotypic evaluations can reveal useful diversity in maize germplasm (Shitta et al., 2021;Chen et al., 2023a).However, morphological characterizations of genotypes are influenced by both the genotype (G) and environment (E), as well as G × E interaction, and hence may not accurately reflect the genotypic effect of nitrogen utilization per se (de Carvalho et al., 2012).In addition, large genotype X season and genotype × location interactions can stymie progress in selecting for low nitrogen tolerance (Ribaut et al., 2007).In this context, understanding the genetic architecture of grain yield (GY) and other secondary traits that are associated with low soil nitrogen can speed up genetic improvement (Ertiro et al., 2020c) for low nitrogen tolerance in tropical maize.
Maize GY under low soil nitrogen is a complex trait governed by multiple genes.In most cases, maize breeding for low nitrogen focuses on the anthesis-silking interval (ASI), anthesis date (AD), plant height (PH), ear position (EPO), and ear height (EH) for indirect selection and GY for direct selection (Worku et al., 2007;Worku et al., 2008;Worku et al., 2012;Das et al., 2019;Ertiro et al., 2020b;Ertiro et al., 2022;Ndlovu et al., 2022).However, morphological characterizations are time-consuming, costineffective, and lack accuracy.The complex nature of breeding maize under low nitrogen conditions has necessitated the use of quantitative trait loci (QTL) based approaches which have been shown to contribute significantly to the understanding of the genetic basis of crop performance and stability under nitrogen-stressed conditions (Semagn et al., 2015).The discovery and characterization of QTL can assist breeders in using genomic regions linked with complex trait expression and deciphering their genetic contribution at the target loci (Ertiro et al., 2020c).In addition, through QTL analysis, the biological mechanisms responsible for phenotypic expression can be pursued.Furthermore, exploiting molecular markers in breeding has allowed the mapping of a subset of markers associated with one or more QTL that contribute to regulating the expression of complex traits like GY.Such molecular markers can form the concrete basis for the use of other genomic approaches such as genomic selection in the desired population (Würschum, 2012).Unfortunately, little is known about the genetic basis of GY and associated traits under low nitrogen conditions, and major effect QTLs are yet to be reported (Ertiro et al., 2020c).
Genomic selection (GS) is a tool increasingly used to reduce the breeding cycle length and increase genetic gain in both crops and livestock (Meuwissen et al., 2001).The success of GS in the dairy industry (where the breeding cycle of the cattle can be reduced from 7 years to 18 months) has helped to achieve twice the genetic gain for key traits (García-Ruiz et al., 2016).Such successes have motivated plant breeders to apply GS in crop improvement for reducing the breeding cycle and increase genetic gain.Currently, GS is applied to improve a range of different traits in maize (Beyene et al., 2015;Ertiro et al., 2020a;Kibe et al., 2020;Gowda et al., 2021;Ndlovu et al., 2022), wheat (Dreisigacker et al., 2021;Bonnett et al., 2022;Juliana et al., 2022;Ficht et al., 2023), and other crops (Hu et al., 2022;Qin et al., 2022;Chen et al., 2023b).GS has been successfully integrated into maize breeding programs for improving GY under optimum and drought conditions (Beyene et al., 2015;Zhang et al., 2017b;Vivek et al., 2017;Beyene et al., 2019;Wang et al., 2020;Beyene et al., 2021).Linkage mapping enables the detection of QTL for the target trait by using different bi-parental populations, whereas GS enables the selection of superior individuals by considering the effects of multiple genes controlling a target trait (Crossa et al., 2017;Yuan et al., 2019).Combining linkage mapping results with GS will accelerate the breeding efficiency for GY and other complex traits under low N stress conditions.Therefore, the full potential of GS needs to be assessed in biparental and/or practical breeding populations for low soil N tolerance.
Studies conducted over the last 20 years have identified QTL that generally explained a significant portion of the phenotypic variance, and therefore gave rise to an optimistic assessment of the prospects of marker-assisted selection (Semagn et al., 2010).According to Coque and Gallais (2006), the detected marker QTL associations in maize revealed the consistency of the involvement of some traits, such as root architecture and glutamine synthase activity, which would be of major importance for GY setting under both optimum and low nitrogen conditions.The discovery of genomic regions associated with GY and other agronomic traits in maize under low soil nitrogen conditions is of paramount importance.Hence, this study sought to evaluate four F3 populations to 1) estimate the phenotypic effect and heritability for GY and other agronomic traits under optimum and low nitrogen management 2) identify the genomic regions associated with these traits, and 3) assessing the usefulness of GS in improving GY and other agronomic traits under optimum and low soil N conditions.The findings of this study can provide genetic resources that can be used in scaling the application of MAS for enhancing maize GY in nitrogen-starved soils in SSA.Furthermore, the identification of QTL can hasten maize breeding cycles and facilitate the release of nitrogen-tolerant or NUE varieties for resource-constrained smallholder farmer communities in the region.

Plant materials
Four F3 tropical maize populations were evaluated in Kenya and Zimbabwe.The populations were developed by the Global Maize Program of the International Maize and Wheat Improvement Center (CIMMYT).The specific details of the tested populations are presented in Table 1.These lines are adapted to mid-altitude regions (1,000-1500 M Above sea level) of SSA.They are bred on a good-by-good basis and adapted to stress conditions.Populations are test crossed with an appropriate tester from the opposite heterotic group for phenotypic evaluations.The CIMMYT lines utilized for the project were derived from breeding programs targeting tolerance to low soil nitrogen; hence the best choice for QTL mapping with most of the lines included in two association mapping panel (AMP) constituted under IMAS (improved maize for Africa soils) and DTMA (Drought Tolerant Maize for Africa) projects (Semagn et al., 2012;Ertiro et al., 2020a).The four populations formed a set of multiple bi-parental populations used in the current study.

Field trial
All four populations were evaluated for response to low soil nitrogen conditions at one to four different locations; Kiboko (Longitude 37 °E, Latitude 2 °S, 975 M A.S.L), Embu (Longitude 37 °E, Latitude 3 °S, 1560 M A.S.L), and Harare in Zimbabwe (Longitude 31 0 E, Latitude 17 °S, 1490 M A.S.L).The lines were evaluated on an alpha lattice incomplete block design under two nitrogen levels.The two nitrogen treatments were low N (Ndepleted field/plot with no application of N fertilizer) and at normal farmer practice conditions (optimum N 200 kg/ha).The low-N sites were prepared by depleting N by growing sorghum at high density with no N fertilizer sources added for four cropping cycles.The depletion crop was uprooted at near maturity and removed from the low-N field to prevent the incorporation of crop residues into the soil.A soil nitrate concentration ranging from 7.5 to 15 parts per million (ppm) is indicative of soil nitrogen deficiency (Zaman-Allah et al., 2018).Low N sites with around 10 ppm of nitrate level provide good experimental conditions for detecting useful genetic variation.The sites used in the present study showed soil nitrate levels ranging from 0.10 to 0.12 ppm in stress experimental sites, while in the optimum sites, the levels were greater than 0.25 ppm.Single row plots measuring 5 m long at 0.75 m row spacing with two seeds per hill were sown.After 3 weeks of planting, plants were thinned to one plant per hill to obtain a final population density of 53,333 plants per hectare.All entries were planted on the same day in conventionally tilled plots and maintained under rainfed or irrigated conditions.

Phenotyping of important agronomic traits
Ten plants in the middle of the row were selected for each genotype for phenotypic evaluation.Phenotypic components measured and analyzed are plant height (PH in centimeters), anthesis date (AD, days), anthesis silking interval (ASI, days), ear height (EH, in centimeters), ear position (EPO, ratio of EH/PH), number of ears per plant (EPP, total number of ears per plot divided by number of plants per plot) and grain yield (GY in t/ha).Mature ears were harvested, manually bagged, air-dried, and shelled on an electric shelling device.The total GY of each plot was weighed on a balance and converted to GY into t/ha.
Best linear unbiased predictors (BLUPs) were calculated with mixed model where genotypes and other factors were treated as random except replications and best linear unbiased estimators (BLUEs) were calculated where genotype entries and replications are treated as fixed effects and the rest of the terms as random.Estimating broad-sense heritability, all the terms were considered random.Broad sense heritability was estimated by the formula: Where σ 2 G is the genotypic variance, σ 2 GE is the Genotypic by environment interaction (GEI), σ 2 e is the error variance, E is the number of environments and r is the number of replications in each trial.The phenotypic and genotypic correlations among traits were evaluated as described by Hirel et al. (2007).

Genotypic analysis
Phenotypic data was collected on testcross hybrids whereas F 3 population lines were genotyped.Since, the single tester was used for each population, assumption is tester effect is same across lines, therefore, tester effect was not included in the model for analyses.Total DNA was extracted from bulked young leaves of the lines according to the CTAB method (CIMMYT, 2005), and the DNA quality for each sample was checked using gel-electrophoresis and spectrophotometer (NanoDrop ND8000 Thermo Scientific) before genotyping.Genotyping was performed using the Illumina MaizeSNP1500 Bead Chip evenly spaced SNP to cover the whole maize genome (Ganal et al., 2011).The above task was performed at LGC genomic labs in the United Kingdom (https://www.lgcgroup.com/genotyping/).Markers which are homozygous for both the parents and polymorphic between them were retained for mapping in each population.Linkage maps in all four populations were constructed using QTL IciMapping version 4.1 software (Meng et al., 2015).Finally, we used 202, 452, 384, and 387 high-quality SNPs in F 3 pop1, F 3 pop2, F 3 pop3, and F 3 pop4, respectively.Linkage map was constructed by using these SNPs, by selecting the most significant markers using stepwise regression.A likelihood ratio test was used to calculate the logarithm of odds (LOD) for each marker at score of >3 with a 30 cM maximum distance between two loci.The recombination frequency between linked loci was transformed into cM (centi Morgan) using Kosambi"s mapping function (Kosambi, 1916).BLUPs across environments were used to detect QTLs based on inclusive composite interval mapping (ICIM) for each population.Phenotypic variation explained by individual QTL and total variation explained by all QTLs together was estimated.QTL naming was done with letter "q" indicating QTL, followed by abbreviation of trait name, the chromosome and marker position, respectively (Ribaut et al., 1997).Additive (a) and dominance (d) effects for each QTL as estimated with QTL IciMapping v.4.1 were used to calculate the ratio of dominance level (|d/a|).This ratio was used to classify the nature of QTL: additive (A; 0 ≤ |d/a| ≤ 0.2); partially dominant (PD; 0.2 < |d/a| ≤ 0.8); dominant (D; 0.8 < |d/a| ≤ 1.2) and overdominant (OD; |d/a| > 1.2).

Genomic prediction
Genomic prediction (GP) analyses were conducted in R program version 4.2.1 (R Core Team 2023).GP was applied on each F 3 population to find out the prediction ability of GY and other agronomic traits evaluated in optimum and low soil N conditions.We used GP model RR-BLUP to predict the untested lines using a five-fold cross validation (Zhao et al., 2012;Crossa et al., 2017).BLUEs across environments for each of the F 3 population were used for the analysis.For GS analyses polymorphic SNPs between the parents of each population was used, like 202, 452, 384, and 387 SNPs in F 3 pop1, F 3 pop2, F 3 pop3 and F 3 pop4, respectively were used.We applied five-fold cross validations with 'within population' approach where both training and estimation set are derived from within each of biparental population.For each trait in each population, 100 iterations were done for sampling of the training and estimation sets.The prediction accuracy was calculated as the correlation between the observed and predicted breeding values divided by the square root of heritability (Dekkers, 2007).

Phenotypic distributions and correlation of traits
GY, ASI, AD, PH, EPO, and EH varied widely across the two nitrogen regimes (i.e., optimum and low nitrogen environments).The extent of variation, however, differed between the four biparental maize populations.Low nitrogen conditions increased the trial mean for AD and decreased the trial mean for GY and PH across the four F 3 populations tested in this study (Figure 1, Supplementary Table S1).GY showed a negative genotypic correlation with AD and ASI across the two nitrogen regimes (Figure 2).Moreover, consistent positive genetic correlations for GY were recorded with EH, EPO, and PH.

Trial mean, genetic variance, and heritability of traits
Population 1 (CML494xCML550): Across the study sites, average GY was 5.22 t/ha and 2.52 t/ha under optimum and low soil nitrogen conditions, respectively (Table 2).A margin of 2-day extension in AD was observed under low soil nitrogen conditions.
Population 2 (CKL5017xCML536): Under optimum conditions (Table 2), moderate heritability was observed for PH (0.68) and EH (0.49) whilst GY (0.28) recorded low heritability across the studied sites.Under low nitrogen conditions, all the studied traits recorded low heritabilities with GY pegged at 0.31.The genotypic and genotype × environment interaction effects were significant (p ≤ 0.05) for GY.Across the studied regions, the mean GY was 7.72 t/ha and 3.37 t/ha under optimum and low nitrogen conditions, respectively.
Population 3 (CML550xCML507): Low heritabilities were recorded across the tested traits under both optimum and low nitrogen conditions (Table 2).GY heritabilities were 0.31 and 0.29 under optimum and low nitrogen conditions, respectively.The genotypic and genotype × environment interaction effects were significant (p ≤ 0.05) for GY.In this population, the mean GY was 6.89 t/ha and 0.94 t/ha under optimum and low nitrogen conditions, respectively.
Four QTLs were identified for ASI under low nitrogen conditions.These were detected on chromosomes 1, 3, 4, and 6.As in the case of GY and AD, no common QTL was identified.TPVE by all QTLs under low nitrogen conditions was 6.91% (population 1) and 16.73% (population 4).
Significant QTLs underlying EPO were recorded as 7 under both nitrogen regimes.One QTL (Chromosome 6: qEPO6_157) was found in both low and optimum nitrogen conditions.Under low soil nitrogen management, QTLs were detected in chromosomes 3, 4, and 6.On the other hand, chromosomes 1, 3, 6, and 8 housed QTL underlying EPO under optimum conditions.Under optimum conditions, the TPVE by all QTL was 15.56%, and 11.13% for populations 1 and 2, respectively.TPVE by all QTLs under low nitrogen conditions was 11.56% (population 1) and 17.98% (population 2).

Overlapping QTL for each trait evaluated under low and optimum nitrogen conditions
The identification of common QTLs is crucial in targeting markers that can be used in breeding for improved nitrogen stress

Genomic prediction correlations for traits evaluated under low and optimum nitrogen conditions
To assess the potential breeding value of GY and other agronomic traits under optimum and low soil N management, GS analysis was performed on all four populations using NUE-associated traits and genotypic markers with RR-BLUP.The RR-BLUP prediction correlations were 0.35, 0.41, 0.06, and 0.35 in populations 1, 2, 3, and 4, respectively (Figure 3; Table 9) under optimum, whereas under low soil N management, the correlations were reduced to 0.21, 0.41, 0.04, and −0.02 in population 1, 2, 3, and 4, respectively.For AD, prediction correlations were 0.31, 0.63, 0.12, and 0.37 under optimum and 0.26, 0.43, 0.29, and 0.36 under low soil N management, respectively (Figure 3; Table 9).For PH, prediction correlations were 0.35, 0.52, 0.34, and 0.18 under optimum, and 0.44, 0.40, 0.15, and 0.10 under low soil N management in populations 1, 2, 3, and 4, respectively.Overall, the prediction correlations were higher in population 2 and lower in population 4 for all the traits under both optimum and low soil N management.

Impact of low nitrogen stress on GY and other associated traits
The role of nitrogen stress in yield reduction and overall crop development cannot be overstated.An in-depth understanding of GY and other related agronomic traits is critical for the evidencebased development of varieties that are tolerant to low soil nitrogen stress (Derera et al., 2008).However, phenotypic characterization of such complex traits is very challenging due to unpredictable   environmental and edaphic conditions (Liu et al., 2010).Importantly, selection through morphological characterization under nitrogen stress conditions is very inconsistent and will not give correct phenotypic data thereby derailing progress in maize breeding schemes.In contrast, the integration of precise genomic tools coupled with conventional breeding can accelerate the development of nitrogen-stress adaptive cultivars with high yields.In this research, two nitrogen regimes (i.e., low, and optimum nitrogen) served as the basis for the identification of QTLs related to GY and other agronomic traits (i.e., AD, ASI, PH, EH, and EPO) and their cross-cutting heritabilities in four F 3 maize populations.The target traits measured followed a normal distribution indicating that the studied traits are quantitative in nature and the tested F 3 populations are stable and suitable for QTL mapping.The mean GY reduction in this study was 69.82% under low nitrogen stress.This is very close to the 71% yield reduction recorded in the study by Ertiro T. B. et al. (2020).Yield reduction under low soil nitrogen stress demonstrates the importance of nitrogen in the growth and development of maize (Ertiro T. B. et al., 2020).Numerous studies have linked a reduction in yield under low-N conditions to decreased kernel number due to abortion (Bänziger et al., 1997;Agrama et al., 1999;Ribaut et al., 2007).Low soil nitrogen environments decrease the number of kernels and ears per plant, decrease the chlorophyll concentration in the ear leaf, and lowers PH by roughly half (Ribaut et al., 2007).
Genotypic and genotype × environment interaction effects were significant for GY across the tested nitrogen regimes.Generally, across the studied environments, the most responsive trait was AD.Our results showed that low nitrogen stress extends the AD in maize.Furthermore, low nitrogen stress increased genetic variance for AD across the studied genotypes.This can be viewed as an indication of the adaptive nature of this secondary trait.
Across the four F 3 populations, there was a consistent moderate to low heritability on all the studied traits across the two nitrogen regimes.For all traits, heritability was shown to decline under low nitrogen stress.The heritabilities of almost all the studied traits were moderate to low.Similar trends were also reported in earlier studies with DH populations (Ertiro et al., 2020b;Ertiro et al., 2022) and association panel (Ertiro et al., 2020a) which were evaluated under optimum and low soil N conditions.Despite this, the heritabilities were significant enough to facilitate indirect selection for increased   Frontiers in Genetics frontiersin.org11 GY under low nitrogen stress.By comparison, PH had low heritabilities under both low and optimum nitrogen conditions.
The pairwise correlation showed a strong negative correlation for GY with AD and ASI across the studied nitrogen regimes.However, other agronomic traits showed positive correlations with GY.These include PH, EH and EPO under both low and optimum soil nitrogen conditions.Similar correlation trends were also reported by earlier studies with maize hybrid and inbred line trials under low soil nitrogen conditions (Ertiro T. B. et al., 2020;Ertiro et al., 2022).PH, EH and EPO had a positive correlation with GY suggesting that improvement of these traits lead to improved maize varieties with high-yielding potentials.Furthermore, these correlations suggest the high odds of tightly linked loci controlling low-N tolerance through the coordinated expression of loci controlling these traits.

QTL mapping under low and optimum nitrogen conditions
NUE traits are highly complex, the advent of high-density marker data has made it feasible to dissect their underlying genetics.The detection of QTLs underlying GY and other associated traits under different nitrogen regimes is integral for scaling breeding initiatives targeting the development of low nitrogen tolerant maize varieties.A total of 91 significant QTLs were identified in this study for six traits under low and optimum nitrogen management (Tables 5, 6, 7, and 8).The disparity in the number of QTLs under optimum (52) and low (39) indicates the effect of genetic variance across the environments.Some of these were found across management conditions and populations.Across the two nitrogen regimes, chromosomes 1 (16) and 8 ( 14) had the highest number of QTLs whilst chromosomes 7 (5) and 5 (6) had the lowest.Chromosomes with a high number of QTLs can be targeted for future research studies focusing on improving GY under low nitrogen stress.Proportional phenotypic variance for the QTLs ranged from 0.7% to 15.22%, with an average of 5.05%.
Population 4 had the highest number of QTLs (36).This is a testament to the diversity within the studied F 3 populations across nitrogen management options.AD (26) and EH (24) had the highest number of QTLs whilst ASI (11) had the lowest.This is consistent with another study (Ertiro et al., 2020c) which found the highest QTL number as those underlying AD.In our study, this was consistent with the observed high genetic variance for AD and EH across the two nitrogen regimes.Marker-assisted selection for GY improvement relies on the successful identification of QTLs with moderate to high effects.For GY, 11 QTLs were detected under both optimum (7) and low (4) nitrogen conditions.This is inconsistent with Ertiro et al., 2020c findings which recorded more QTLs under low nitrogen management than under optimum conditions.One  The distinct genetic controls of the expression of the phenotypes reported in this study and related studies (Bänziger et al., 1999;Ribaut et al., 2007) suggest that selection for low nitrogen stress tolerance will be more effective under low nitrogen management.These findings unravel the most promising genomic regions for MAS to increase maize tolerance to low nitrogen stress.Further multi-environment trials with larger population sizes are necessary to validate the stability of the recorded QTLs and utilize them to identify underlying causal genes.
To achieve higher genetic gain, complex traits in maize breeding can be improved by integrating genomic, bioinformatic, and statistical tools into breeding programs (Beyene et al., 2019;Beyene et al., 2021).Rapid technological advancements have made available cheaper genotyping tools, while rapid progress in the field of big data and bioinformatics has led to the development of user-friendly software tools which can handle complex statistical models and facilitates breeders to apply them to improve complex traits on a more routine basis.GS, which predicts the breeding values by employing genome-wide markers is proving to be effective in improving complex traits that are controlled by multiple minor effect genes.For instance, RR-BLUP or G-BLUP is now very commonly used to predict several complex traits in maize (Guo et al., 2012;Hao et al., 2019;Beyene et al., 2021).Several NUE traits like chlorophyll index, chlorophyll fluorescence and leaf N content were also selected by genomic prediction before integrating into NUE breeding programs in ryegrass (Zhao et al., 2020).In the present study, GS was performed on GY and other agronomic traits evaluated under optimum and low soil N stress conditions.The average genomic prediction correlations ranged from 0.07 to 0.41 under optimum and −0.07 to 0.41 under low soil N stress conditions.The higher values of prediction correlations are comparable to earlier studies with a diversity panel evaluated under optimum and low soil N stress conditions (Ertiro et al., 2020c).The observed prediction accuracies for other agronomic traits are comparable to earlier studies reported under different stresses in maize (Zhang et al., 2015;Yuan et al., 2019).
In population 3, we found low accuracy for GY under optimum condition, which could be due to its small range of variability within the population as well as low heritability (Table 2, Figure 1).We observed a wide range of prediction correlations for same trait among different populations and management.This could be due to their differences in sample size, genetic variance, trait heritability, changes in population structure and linkage disequilibrium estimates.In some populations, we observed negative prediction correlations like GY for pop 4 under low soil N stress conditions (Figure 3; Table 9).Opposite linkage phases between markers and major-effect QTLs in the population may be the reason for negative correlations.In addition, the four F 3 populations used in this study were developed by using tolerant x tolerant crosses for low soil N stress, which may be another reason why we observed low prediction correlations (as most of the causal factors for these traits might have been fixed on both the parental lines).As a result, variation is low for the traits which is evident with low to moderate heritability estimates (Table 2).Nevertheless, the prediction correlation followed a consistent pattern for all traits under both optimum and low N conditions.Traits with high prediction correlations also tend to have relatively high heritability estimates.In GS, the less complex trait AD and PH had higher accuracy compared to GY, which is consistent with the nature of trait complexity (Zhang et al., 2017a;Yuan et al., 2019).For simple inherited traits that are positively correlated with GY, prediction correlations are moderate to high, which clearly supports the usefulness of GS for their improvement under either optimum or low soil N stress conditions.
Overall, the linkage mapping studies with our F 3 populations revealed several QTL with minor to moderate effects for GY and other agronomic traits under optimum and low N management.It is difficult to capture multiple QTL with small to moderate effects for selection and their environmental and genotypic specific expressions make it even more difficult to improve these traits only through traditional breeding or a few QTL based MAS.However, discovery of genomic regions through linkage mapping will continue to be vital to understanding the genetic basis of these traits.On the other hand, genome-wide selection is critical in improving quantitative traits.The phenotypic selection efficiency per cycle is measured as h (square root of heritability) and the value of marker-based prediction accuracy close to h indicates the selection response based on markers and based on phenotypes are near equal (Dekkers, 2007;Lorenzana and Bernardo, 2009).The prediction accuracies observed for GY and other traits was ≥1/2h indicating that the response to marker-based selection would be at least 50% of the response to phenotypic selection for agronomic traits.With the possibility of making three cycles per year by using marker-based selection, the selection response will be −1.5 times the gain from one cycle of phenotypic selection.These results indicate that genomewide selection would be more efficient in terms of genetic gain per year.However, one must be cautious as we observed prediction accuracy of <0.10 for some traits in some populations, where GS has no additional advantage over phenotypic selection.Frontiers in Genetics frontiersin.org15

Conclusion
Nitrogen-depleted soils are a major factor behind low maize productivity in smallholder farming systems in SSA.Genomicsbased plant breeding techniques such as QTL mapping and GS can provide useful information for scaling the development of low nitrogen tolerant maize varieties.Here, we sought to identify the genomic regions associated with GY and related traits under optimum and low soil nitrogen conditions in F 3 maize populations grown in Kenya and Zimbabwe.Our analysis found a total of 91 QTLs underlying GY, ASI, AD, EH, EPO, and PH.However, no common QTL for GY was identified across the studied nitrogen regimes.The validation of these QTLs in the same F 3 populations is needed to guarantee the success of future markerassisted selections.Identification of many QTL with minor effects indicates rather QTL mapping on these traits, GS is more appropriate for their improvement.Genomic prediction correlations were low to moderate.However, by considering the possibility to have three cycles per year with marker-based selection, integration of GS in a low N tolerant breeding program will be more efficient to improve GY under low soil N conditions.

FIGURE 2
FIGURE 2Phenotypic correlations of GY and other agronomic traits evaluated under optimum and low soil nitrogen conditions in F 3 pop1 CML494 x CML550 population.The x marks indicated values are not significant at p < 0.05.LN-low soil N conditions.Correlations with >0.11, and >0.19 were significant at 0.05, and 0.01 (p) levels, respectively.GY-Grain yield; AD -anthesis date; ASI -anthesis silking interval; PH-plant height; EH, ear height; Epos-ear position; EPP-ears per plant; SEN-senescence.
QTL for GY (qGY10_130) was found under optimum conditions on chromosome 10 in both population.A study byErtiro et al.,  2020c  also found QTLs underlying GY under different nitrogen regimes on chromosome 10.No common QTL for GY was identified across the studied F 3 populations.Ertiro et al. (2020a) also observed similar results, i.e., no common QTLs observed for GY.Noteworthy, the use of a different number of populations, locations and markers in future studies can give rise to a different outcome.In this study, some common QTLs were identified for AD, EH and EPO.The identification of common QTLs justifies the magnitude of correlation across the two nitrogen regimes.Common QTLs can also be used for indirect selection.Similar to the results ofErtiro et al., 2020c, QTL correspondences detected between the nitrogen regimes for the studied traits were not similar across F 3 populations and this points to the genetic-background-specific nature of these QTLs.

FIGURE 3
FIGURE 3Distribution of the five-fold cross-validated genomic prediction correlations for grain yield (GY), anthesis date (AD), anthesis-silking interval (ASI), plant height (PH), ear height (EH), and ear position (EPO) evaluated in multiple environments under optimum and low soil N management conditions in four F 3 populations.

TABLE 1
Details of the populations used in this study and number of locations evaluated under optimum and low N management.

TABLE 2
Trait means, heritability and variance components of GY and other traits evaluated under optimum and low nitrogen conditions in four F 3 populations across multiple environments.

TABLE 2 (
Continued) Trait means, heritability and variance components of GY and other traits evaluated under optimum and low nitrogen conditions in four F 3 populations across multiple environments.andNUE.This study identified several common QTLs across the studied F3 populations and nitrogen regimes (Tables4-8).Under optimum conditions, one QTL (qAD03_158) for AD was overlapping with one QTL under low soil nitrogen conditions on chromosome 3 (157.97-159.08Mbp).For ASI and PH, no common QTL was identified across the two nitrogen regimes. tolerance

TABLE 3
Linkage map information for the four F3 populations used for QTL analysis.

TABLE 4
Genetic characteristics of QTLs detected for GY and associated traits under optimum and low-nitrogen stress in F 3 population 1.
Opt, Optimum; Low N-low soil N stress management; GY, Grain yield; AD, anthesis date; ASI, anthesis silking interval; PH, plant height; EH, ear height; EPO, ear position.The italic values refer to the names of identified QTLs.

TABLE 5
Genetic characteristics of detected QTLs for GY and associated traits under optimum and low-nitrogen stress in F3 population 2.

TABLE 5 (
Continued) Genetic characteristics of detected QTLs for GY and associated traits under optimum and low-nitrogen stress in F3 population 2. Low N-low soil N stress management; GY, Grain yield; AD, anthesis date; ASI, anthesis silking interval; PH, plant height; EH, ear height; EPO, ear position.The italic values refer to the names of identified QTLs.

TABLE 6
Genetic characteristics of detected QTLs for GY and associated traits under optimum and low-nitrogen stress in F3 population 3.
Opt, Optimum; Low N-low soil N stress management; GY, Grain yield; AD, anthesis date; PH, plant height; EH, ear height.The italic values refer to the names of identified QTLs.

TABLE 7
Genetic characteristics of detected QTLs for GY and associated traits under optimum and low-nitrogen stress in F3 population 4.

TABLE 8
Summary of detected QTL for measured traits under optimum and low-nitrogen stress in F3 populations derived from seven elite lines.

TABLE 9
The prediction accuracy for grain yield and other traits evaluated under optimum and low-nitrogen stress conditions for four F 3 populations derived from seven elite lines.Low N-low soil N stress management; GY, Grain yield; AD, anthesis date; ASI, anthesis silking interval; PH, plant height; EH, ear height; EPO, ear position.